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We analyze the phenomenon of preheating, i.e. explosive particle produc- 
^vq ■ tion due to parametric amplification of quantum fluctuations in the unbroken 

symmetry case, or spinodal instabilities in the broken symmetry phase, using 
\^ | the Minkowski space 0(N) vector model in the large N limit to study the 

0^ ' non-perturbative issues involved. We give analytic results for weak couplings 

and times short compared to the time at which the fluctuations become of the 
same order as the tree level terms, as well as numerical results including the 
Oh. full backreaction. In the case where the symmetry is unbroken, the analyti- 

cs ! cal results agree spectacularly well with the numerical ones in their common 

domain of validity. In the broken symmetry case, interesting situations, cor- 
J^j ■ responding to slow roll initial conditions from the unstable minimum at the 

origin, give rise to a new and unexpected phenomenon: the dynamical relax- 
ation of the vacuum energy. That is, particles are abundantly produced at the 
expense of the quantum vacuum energy while the zero mode comes back to 
almost its initial value. In both cases we obtain analytically and numerically 
the equation of state which in both cases can be written in terms of an ef- 
fective polytropic index that interpolates between vacuum and radiation-like 
domination. 

We find that simplified analysis based on harmonic behavior of the zero 
mode, giving rise to a Mathieu equation for the non-zero modes miss im- 
portant physics. Furthermore, such analysis that do not include the full 
backreaction and do not conserve energy, results in unbound particle pro- 
duction. Our results rule out the possibility of symmetry restoration by non- 
equilibrium fluctuations in the cases relevant for new inflationnary scenarios. 
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Finally, estimates of the reheating temperature are given, as well as a dis- 
cussion of the inconsistency of a kinetic approach to thermalization when a 
non-perturbatively large number of particles are created. 

11.10.-z;11.15.Pg;98.80.Cq 
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I. INTRODUCTION 



It has recently been realized that as the zero momentum mode of a quantum field 
evolves, it can drive a large amplification of quantum fluctuations. This, in turn, gives rise 
to copious particle production for bosonic fields, creating quanta in a highly non-equilibrium 
distribution, radically changing the standard picture of reheating the post-inflationary uni- 
verse This process has other possible applications, such as in understanding the 
hadronization stage of the quark-gluon plasma J7| as well as trying to understand out of 
equilibrium particle production in strong electromagnetic fields and in heavy ion collisions 

The actual processes giving rise to preheating can be different depending on the potential 
for the scalar field involved as well as the initial conditions. For example, in new inflationary 
scenarios, where the inflaton field's zero mode evolves down the flat portion of a potential 
admitting spontaneous symmetry breaking, particle production occurs due to the existence 
of unstable field modes which get amplified until the zero mode leaves the instability region. 
These are the instabilities that give rise to spinodal decomposition and phase separation. 
In contrast, if we start with chaotic initial conditions, so that the field has large initial am- 
plitude, particles are created from the parametric amplification of the quantum fluctuations 
due to the oscillations of the zero mode. 

In this paper we analyze the details of this so-called preheating process both analytically 
as well as numerically. Preheating is a non-perturbative process, with typically 1/A particles 
being produced, where A is the self coupling of the field. Due to this fact, any attempts at 
analyzing the detailed dynamics of preheating must also be non-perturbative in nature. This 
leads us to consider the O(N) vector model in the large N limit. This is a non-perturbative 
approximation that has many important features that justify its use: unlike the Hartree or 
mean-field approximation ||, it can be systematically improved in the 1/N expansion. It 
conserves energy, satisfies the Ward identities of the underlying symmetry, and again unlike 
the Hartree approximation it predicts the correct order of the transition in equilibrium. 

This approximation has also been used in other non-equilibrium contexts In this 

work, we consider this model in Minkowski space, saving the discussion of the effects of the 
expansion of the universe for later work. 

Our findings are summarized as follows. 

We are able to provide consistent non-perturbative analytic estimates of the non- 
equilibrium processes occurring during the preheating stage taking into account the exact 
evolution of the inflaton zero mode for large amplitudes when the quantum back-reaction 
due to the produced particles is negligible i.e. at early and intermediate times. We also 
compute the momentum space distribution of the particle number as well as the effective 
equation of state during this stage. Explicit expressions for the growth of quantum fluctu- 
ations, the preheating time scale, as well as the effective (time dependent) polytropic index 
defining the equation of state are given in sec. III. 

We then go beyond the early /intermediate time regime and evolve the equations of mo- 
tion numerically, taking into account back-reaction effects. (That is, the non-linear quantum 
field interaction). These results confirm the analytic results in their domain of validity and 
show how, when back-reaction effects are large enough to compete with tree level effects, 
dissipational effects arise in the zero mode. Energy conservation is guaranteed in the full 
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backreaction problem, leading to the eventual shut-off of particle production. This is an 
important ingredient in the dynamics that determines the relevant time scales. 

We also find a novel dynamical relaxation of the vacuum energy in this regime when 
the theory is in the broken phase. Namely, particles are produced at the expense of the 
quantum vacuum energy while the zero mode contributes very little. We find a radiation 
type equation of state for late times {p ~ 3 e) despite the lack of thermal equilibrium. 

Finally we discuss the calculation of the reheating temperature in a class of models, 
paying particular attention to when the kinetic approach to thermalization and equilibration 
is applicable. 

There have been a number of papers (see refs. |l].[2|.[T2j - dedicated to the analysis of 



the preheating process where particle production and back-reaction are estimated in different 



approximations 29 



The layout of the paper is as follows. Section II presents the model, the evolution equa- 
tions, the renormalization of the equations of motion and introduces the relevant definitions 
of particle number, energy and pressure and the details of their renormalization. The unbro- 
ken and broken symmetry cases are presented in detail and the differences in their treatment 
are clearly explained. 

In sections III through V we present a detailed analytic and numerical treatment of 
both the unbroken and broken symmetry phases emphasizing the description of particle 
production, energy, pressure and the equation of state. In the broken symmetry case, when 
the inflaton zero mode begins very close to the top of the potential, we find that there is 
a novel phenomenon of relaxation of the vacuum energy that explicitly shows where the 
energy used to produced the particles comes from. We also discuss why the phenomenon of 
symmetry restoration at preheating, discussed by various authors 0,[13],|25],|26| is not seen to 



occur in the cases treated by us here and relevant for new inflationnary scenarios 



In section VI we provide estimates, under suitably specified assumptions, of the reheating 
temperature in the O(N) model as well as other models in which the inflaton couples to 
lighter scalars. In this section we argue that thermalization cannot be studied with a kinetic 
approach because of the non-perturbatively large occupation number of long- wavelength 
modes. 

Finally, we summarize our results and discuss future avenues of study in the conclusions. 
We also include an appendix where we gather many important technical details on the 
evaluation of the Floquet mode functions and Floquet indices used in the main text. 



II. SCALAR FIELD DYNAMICS IN THE LARGE N LIMIT 

As mentioned above, preheating is a non-perturbative phenomenon so that a non- 
perturbative treatment of the field theory is necessary. This leads us to consider the O(N) 
vector model in the large N limit. 

In this section we introduce this model, obtain the non-equilibrium evolution equations, 
the energy momentum tensor and analyze the issue of renormalization. We will then be 
poised to study each particular case in detail in the later sections. 

The Lagrangian density is the following: 

C=^d li $-& t $-V($-$) 
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V(a,Tt) = -m 2 f 0+^(0 4? (2.1) 

for A fixed in the large N limit. Here is an 0(N) vector, <fi = (a, tt) and n represents the 
N — 1 "pions". In what follows, we will consider two different cases of the potential V(a, tt), 
with (m 2 < 0) or without (m 2 > 0) symmetry breaking. 

We can decompose the field a into its zero mode and fluctuations x(x,t) about the zero 
mode: 

a(x,t) = a (t)+ X (x,t) (2.2) 

The generating functional of real time non-equilibrium Green's functions can be written 
in terms of a path integral along a complex contour in time, corresponding to forward and 
backward time evolution and at finite temperature a branch down the imaginary time axis. 
This requires doubling the number of fields which now carry a label ± corresponding to 
forward (+), and backward (— ) time evolution. The reader is referred to the literature 
for more details fpL4|,pL5l| - This generating functional along the complex contour requires the 
Lagrangian density along the contour, which for zero temperature is given by || 




The tadpole condition (x ± (x,t)) = will lead to the equations of motion as discussed in || 
and references therein. 

A consistent and elegant version of the large N limit for non-equilibrium problems can 
be obtained by introducing an auxiliary field and is presented very thoroughly in reference 
||. This formulation has the advantage that it can incorporate the 0(1/ N) corrections in a 
systematic fashion. Alternatively, the large iV limit can be implemented via a Hartree-like 
factorization [[| in which i) there are no cross correlations between the pions and sigma field 
and ii) the two point correlation functions of the pion field are diagonal in the 0(N—1) space 
of the remaining unbroken symmetry group. To leading order in large N both methods are 
completely equivalent and for simplicity of presentation we chose the factorization method. 

The factorization of the non-linear terms in the Lagrangian is (again for both ± compo- 
nents): 



X 4 6(x 2 )x 2 + constant (2.4) 

X 3 - 3( X 2 )X (2.5) 

(vf ■ tt) 2 2(7r 2 )7r 2 - (tt 2 ) 2 + 0(1/N) (2.6) 

vfV-^V + ^X 2 ) (2.7) 

t\ - (v? 2 )x (2.8) 



To obtain a large N limit, we define 
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N-l 

Tf(x,t) = ■•■,!) ; <ro(t) = c(>(t)y/N 

where the large N limit is implemented by the requirement that 



(2.9) 



(2.10) 



The leading contribution is obtained by neglecting the 0(1/ N) terms in the formal large N 
limit. The resulting Lagrangian density is quadratic, with a linear term in \ '■ 



£[o"o + X + , - £K + X ,7? ] = I 



\{d,x + ) 2 + \{d^) 2 -x + V(t) 



-\Ml(t)(x + f- l -Ml(t)(r 



{( 



X 



X 



7T 



7T 



-)} 



where, 



V(<f>(t),t) = VN<f>(t) 



A 



m 2 + ^0 2 (t) + ^ 2 (t)> 



A 
2* 
A 



^ 2 (t) = m 2 + -0 2 (t) + -^ 2 (t)> 



3A 



A 



M{(t) = m' + -<P\t) + -(^(t)). 

Note that we have used spatial translational invariance to write 

(^(x,t)) = (^(t)) 



(2.11) 

(2.12) 

(2.13) 
(2.14) 

(2.15) 



The necessary (zero temperature ) non-equilibrium Green's functions are constructed from 
the following ingredients 



G>(t,t') 
G$(t,t>) 



2W k 

i 

2Wk 



Vk(t)V k *(t') 
V k (t')V k *(t) 



(2.16) 
(2.17) 



while the Heisenberg field operator if)(x, t) can be written as 

'a k V k (ty % - s + a\v:(t)e- i% ^ 



ip(x,t) = -=Vl 



(2.18) 



where a k , a\ are the canonical destruction and annihilation operators and V the quantization 
volume. 

The evolution equations for the expectation value <p(t) and the mode functions V k (t) can 
be obtained by using the tadpole method B and are given by: 
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4>(t) + <j>{t) 



m 2 + ^ 2 (t) + ^(t)} 



^ 2 +k 2 + rn 2 + \m + \mt)) 



V k (t) = ; V k (0) = 1 ; 14(0) 



W k = Jk 2 + ml 



d 3 k \V k (t)f 



(2tt) 3 2W k 



(2.19) 

-iW k (2.20) 

(2.21) 
(2.22) 

The initial state \i) is chosen to be the vacuum for these modes, i.e. a k \i) = 0. The 
frequencies W k (i.e. m§) will determine the initial state and will be discussed for each 
particular case below. 

The fluctuations x(x, t) obey an independent equation, that does not enter in the dy- 
namics of the evolution of the expectation value or the 7? fields to this order and decouples 
in the leading order in the large N limit [§. 

It is clear from the above equations that the Ward identities of Goldstone's theorem are 
fulfilled. Because V' (<j>(t),t) = yfN<j)(t)M%{t), whenever V'(<j>(t),t) vanishes for ^ then 
Aijf = and the "pions" are the Goldstone bosons. This observation will be important in 
the discussions of symmetry breaking in a later section. 

Since in this approximation, the dynamics for the 7? and x fields decouple, and the 
dynamics of \ does not influence that of 0, the mode functions or (ip 2 ), we will only con- 
centrate on the solution for the 7? fields. We note however, that if the dynamics is such that 
the asymptotic value of <p 7^ the masses for x and the "pion" multiplet 7? are different, and 
the original O(N) symmetry is broken down to the 0(N — 1) subgroup. 



A. Renormalization of the O(N) Model 

We briefly review the most relevant features of the renormalization program in the large 
N limit that will be used frequently in our analysis. For more details the reader is referred 
to 

In this approximation, the Lagrangian is quadratic, and there are no counterterms. This 
implies that the equations for the mode functions must be finite. This requires that 

m 2 + ^(t) + ^ 2 (t)> = m\ + ^0 2 (t) + ^\t)) R = -v(t) (2.23) 



Defining 



U</) -ip k (t) , <p k {0) = -±= , ^ k {0) = -iJw k (2.24) 



VW k rns ' r ^ ' vw k 
The function tp k (t) is written as linear combinations of WKB solutions of the form 

cp k (t) = A k exp f R k (t')dt' + B k exp f* Rt(t')dt' (2.25) 
Jo Jo 

with R k (t) obeying a Riccati equation [[LlJ and the coefficients A k , B k are fixed by the initial 
conditions. After some algebra we find 
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<pk(t) 



2 fc^oo 1 V(t) 3v(t) 2 -v(t) _ 7 v 

1 " k + 2~¥ + + ° {k } 



oscillatory terms 



12 fc -°° k - + ^\ ,! (t)2 + 0(AT 5 ) + oscillatory terms 



2k ' 8k 3 

Using this asymptotic forms, we obtain PJlTj the following renormalized quantities 



2.26) 



m 2 + 



A 



167T 2 

Ah 
16tt 2 



A 2 - rrz 2 D In f -) 



In 



A 



A« 



r/7 



R 



and 



R 



k 2 dk 



1^(01-7 + 



i e(jfe-«) 



A 



A; 



R 



(2.27) 
(2.28) 

(2.29) 



k ' 2k 3 

where we have introduced the (arbitrary) renormalization scale k. Eqs. [2.23| and |2.29| lead 
to the renormalization conditions valid in the large N limit. 

At this point it is convenient to absorb a further finite renormalization in the definition 
of the mass and introduce the following quantities: 



A 



M 2 R = ml+'f(i> 2 (0)) R 



r = \M R \t , q 



M. 



w. 



R\ 



M 



R\ 



2\M R \ 
Ar 
" 2\Mr\ 2 



mt)) R - (^ 2 (o)) 



R 



(S(0)=0) 



<p q (r) = \M R \ (p k (t) 



(2.30) 
(2.31) 

(2.32) 

(2.33) 

(2.34) 
(2.35) 



For simplicity in our numerical calculations later, we will chose the renormalization scale 
K = \M R \. The evolution equations are now written in terms of these dimensionless variables, 
in which dots now stand for derivatives with respect to r. 



B. Unbroken Symmetry 

In this case M R = \M R \ 2 , and in terms of the dimensionless variables introduced above 
we find the following equations of motion: 

fi + ri + ri 3 + gri(r)^(r) = (2.36) 
d 2 



dr 2 



+ q 2 + l+ r ? (r) 2 + ( 7 S(r) 
1 



<Pq{ 



T 



o, 



ip q (0) = -% 



7/(0) =770 , 7)(0) =0 



(2.37) 
(2.38) 
(2.39) 
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As mentioned above, the choice of Q, q determines the initial state. We will choose these 
such that at t = the quantum fluctuations are in the ground state of the oscillators at the 
initial time. Recalling that by definition g£(0) = 0, we choose the dimensionless frequencies 
to be 



The Wronskian of two solutions of (|2.37|) is given by 

W [(p q , <p q ] = 2i , 

while (?S(r) is given by 



r) 



1 



O(q-l) 
2g 3 



-Va+V 2 {T) + gX(r) 



[2A0) 



(2.41) 



^2.42) 



In the case of broken symmetry M\ 
limit become: 



T) — 7] + rf + g T](t) £(t) = 
d 



C. Broken Symmetry 

— | Mr | and the field equations in the N = oo 



— + q 2 -l+ rj(ry + gE(T) 



dr 2 



(2.43) 
(2.44) 



where S(r) is given in terms of the mode functions v 9 g( r ) by the same expression of the 
previous case, (|2.42|). Now the choice of boundary conditions is more subtle. The situation 



of interest is when < rfc < 1, corresponding to the situation where the expectation value 
rolls down the potential hill from the origin. The modes with q 2 < 1 — i]q are unstable and 
thus do not represent simple harmonic oscillator quantum states. Therefore one must chose a 
different set of boundary conditions for these modes. Our choice will be that corresponding 
to the ground state of an upright harmonic oscillator. This particular initial condition 
corresponds to a quench type of situation in which the initial state is evolved in time in an 
inverted parabolic potential (for early times t > 0). Thus we shall use the following initial 
conditions for the mode functions: 



1 



(p g (0) = ~i Jttq 



l+r]$ for q < 1 - r] 



= Jq 2 - 1 + Vo for g 2 > 1 



V 2 o 



< Vo < 1 



(2.45) 

(2.46) 
(2.47) 



along with the initial conditions for the zero mode given by eq. (|2.39| ). 
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D. Particle Number 



Although the notion of particle number is ambiguous in a time dependent non-equilibrium 
situation, a suitable definition can be given with respect to some particular pointer state. 
We consider two particular definitions that are physically motivated and relevant as we will 
see later. The first is when we define particles with respect to the initial Fock vacuum state, 
while the second corresponds to defining particles with respect to the adiabatic vacuum 
state. 

In the former case we write the spatial Fourier transform of the fluctuating field ip(x,t) 
in (|2.9|) and its canonical momentum U(x,t) as 



(2.48) 



[2.49) 



with the time independent creation and annihilation operators, such that a k annihilates the 
initial Fock vacuum state. Using the initial conditions on the mode functions, the Heisenberg 
field operators are written as 



u k (t) = u-\t) n fc (o)w(t) 

a k (t) =U-\t) a k U{t) 




a k (t)+al k (t) 
ak(t) - aL k (t) 



(2.50) 

(2.51) 
(2.52) 



with U(t) the time evolution operator with the boundary condition U(0) = 1. The Heisen- 
berg operators atit) , a\(t) are related to a k , a\ by a Bogoliubov (canonical) transformation 
(see reference M for details). 

The particle number defined with respect to the initial Fock vacuum state is defined in 
term of the dimensionless variables introduced above as 



N g (r) 



(al(t)a k (t)) = 



a, 



w + 



Q 



1 

2' 



^2.53) 



It is this definition of particle number that will be used for the numerical study. 

In order to define the particle number with respect to the adiabatic vacuum state we note 
that the mode equations fl2.37|J2.44]) are those of harmonic oscillators with time dependent 
squared frequencies 



uj 2 Jr) = q 2 ±l + r ] 2 (T)+g^(r) 



(2.54) 



with + for the unbroken symmetry case and — for the broken symmetry case, respectively. 
When the frequencies are real, the adiabatic modes can be introduced in the following 
manner: 
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n fc (t) 



2u k {t) 



lLO k {t) 



[2.55) 
[2.56) 



where now a k (t) is a canonical operator that destroys the adiabatic vacuum state, and is 
related to a k , a\ by a Bogoliubov transformation. This expansion diagonalizes the instan- 
taneous Hamiltonian in terms of the canonical operators a(t) , a^(t). The adiabatic particle 
number is 



Mt)| 2 + 



m\T)\' 



(2.57) 



As mentioned above, the adiabatic particle number can only be defined when the fre- 
quencies tOkit) are real. Thus, in the broken symmetry state they can only be defined 
for wave- vectors larger than the maximum unstable wave- vector, k > k u = |Mr| 



These adiabatic modes and the corresponding adiabatic particle number have been used 
previously within the non-equilibrium context [|8"Hl0f and will be very useful in the analysis 
of the energy below. Both definitions coincide at t — because u) q (fi) = Q q . Notice that 
Nf(0) = N q {0) = due to the fact that we are choosing zero initial temperature. (We 



considered a non-zero initial temperature in refs. P,[TT 



E. Energy and Pressure 

The energy momentum tensor for this theory is given by 



= d*<j)-d l <p-g 
Using the large N factorization 



d a (j> ■ d a <p - V{(p 



(2.58) 



E 

NV 



+ linear terms in ip + 0(1/N) 



X 



X 



u;i(t) = e + m 2 + - ( j ) 2 (t) + -(r(t)}. 



2~ty ) we find the energy density operator to be 

(2.59) 
(2.60) 



Taking the expectation value in the initial state and the infinite volume limit, defining 
e =< E > /NV, and recalling that the tadpole condition requires that the expectation 
value of ip vanishes, we find the expectation value of the energy to be 



e = -<p 2 (t) + -m 2 2 (t) + -0 4 (t) + J- 
2^ w 2 r w 8^ W 8tt 2 



k 2 dk\\Mt)\ 2 + uj 2 k (t)\Mt)\ 2 } -^ 2 W> 2 (2.61) 



It is now straightforward to prove that this bare energy is conserved using the equations 



of motion (|2.19| )- (p.21|) . It is important to account for the last term when taking the time 
derivative because this term cancels a similar term in the time derivative of oJ k {t). 



11 



Since we consider translationally as well as rotationally invariant states, the expectation 
value of T^ v takes the fluid form 



T 



oo 



-< E > 



=<P> 8 13 



->i0 







^2.62) 



We want to emphasize that the full evolution of the zero mode plus the back-reaction 
with quantum fluctuations conserves energy. Such is obviously not the case in most treat- 
ments of reheating in the literature in which back-reaction effects on the zero mode are 
neglected. Without energy conservation, the quantum fluctuations grow without bound. In 
cosmological scenarios energy is not conserved but its time dependence is not arbitrary; in 
a fixed space-time background metric it is determined by the covariant conservation of the 
energy momentum tensor. There again only a full account of the quantum back-reaction 
will maintain covariant conservation of the energy momentum tensor. 

We can write the integral in eq. ( |2.61| ) as 



1 



8tt 2 



k 2 dk 



\Mt)\ 2 +ul(t)\<Pk(t)f 



eu + 



2 7T 2 h u k2dkU ° k{t) ( N ^ + \ 



■S7T 2 JO 



k 2 dk 



\^ k {t)\ 2 + u 2 k {t) \<p k {t)\' 



(2.63) 
(2.64) 



where A is a spatial upper momentum cutoff, taken to infinity after renormalization. In 
the broken symmetry case, eu is the contribution to the energy-momentum tensor from the 
unstable modes with negative squared frequencies, k 2 = |Mr| 2 [1 — t]q] and N% d (t) is the 
adiabatic particle number given by eq. (|2.57|) . For the unbroken symmetry case Ejj = and 
k u = 0. 

This representation is particularly useful in dealing with renormalization of the energy. 
Since the energy is conserved, a subtraction at t — suffices to render it finite in terms 
of the renormalized coupling and mass. Using energy conservation and the renormalization 
conditions in the large N limit, we find that the contribution k 2 dk U0k{t) N% d (t) is finite. 
This also follows from the asymptotic behaviors ( |2.26|) . 

In terms of dimensionless quantities, the renormalized energy density is, after taking 
A -> oo: 



2|Mp| 4 



A 



R 



V j + V jM 2 (r)+s F 




A; 



\M R \ 2 



< ^ 2 (0) > R M 2 {t) 



+ 
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M\ t) 
4 



+ M 2 (r)q u ^ql + M 2 (r) - 2q u (q 2 u + M 2 (r 



3/2 



+ M\t) In 

e F = ~ / q dq 
2 Jo 



q u + Jq 2 u + M 2 (r) 



\<p g \ 2 +UJ 2 (T)\<f q \ 2 



+ c 

rod 

+ 2g q 2 dqu q {r)N a q d {r 

Jq-u, 



M 2 {r) = ±l + r j 2 {r)+gi:{r) , uj 2 {t) = q 2 + M 2 (r) 



(2.65) 

(2.66) 
(2.67) 



where the lower sign and q u = J 1 — i] 2 apply to the broken symmetry case while the upper 
sign and q u = correspond to the unbroken symmetry case. The constant C is chosen such 
that e coincides with the classical energy for the zero mode. The quantity yVt(r) is identified 
as the effective (dimensionless) mass for the "pions". 
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We find using the renormalized eqs. (gjg ), Q), ( ^42|) , (^43|) and (ggg) , that the 
renormalized energy e is indeed conserved both for unbroken and for broken symmetry. 

The pressure is obtained from the spatial components of the energy momentum tensor 
(see eq. (|2.62p ) and we find the expectation value of the pressure density p =< P > /NV to 
be given by 



V 



+ 



4tt 2 



k dk 



|^(t)| 2 + yK(t)| 2 



— e 



[2.68) 



Using the large-/c behavior of the mode functions (|2.26| ), we find that aside from the time 
independent divergence that is present also in the energy the pressure needs an extra sub- 
traction |0 2 //c 3 compared with the energy. Such a term corresponds to an additive renor- 
malization of the energy-momentum tensor of the form 



(2.69) 



with A a (divergent) constant [16]. Performing the integrals with a spatial ultraviolet cutoff, 
and in terms of the renormalization scale k introduced before, we find 



A = -—hx[-] 
12 L fc J 



(2.70) 



In terms of dimensionless quantities and after subtracting a time independent quartic 
divergence, we finally find setting A = oo, 



2|Mr 14 



P 



+ 



\f] 2 + g f q 2 dq 
X R { Jo 

%~1) d 2 , Jl 



jMr)\ 2 + \^(r)\ 2 -^q-^-M 2 (r) 



e . 



(2.71) 



At this stage we can recognize why the effective potential is an irrelevant quantity to 
study the dynamics. 

The sum of terms without Ef in ( 2.65 ) for q u = are identified with the effective potential 
in this approximation for a time independent r] ; gH. These arise from the "zero point" energy 



of the oscillators with time dependent frequency in (|2.63| ). 

In the broken symmetry case the term ep describes the dynamics of the spinodal in- 



stabilities [[Tl[] since the mode functions will grow in time. Ignoring these instabilities and 
setting q u = as is done in a calculation of the effective potential results in an imaginary 
part. In the unbroken symmetry (q u = 0) case the sum of terms without Ep give the effective 
potential in the large N limit, but the term ep describes the profuse particle production via 
parametric amplification, the mode functions in the unstable bands give a contribution to 
this term that eventually becomes non-perturbatively large and comparable to the tree level 
terms as will be described in detail below. Clearly both in the broken and unbroken symme- 
try cases the effective potential misses all of the interesting non-perturbative dynamics, that 
is the exponential growth of quantum fluctuations and the ensuing particle production, ei- 
ther associated with unstable bands in the unbroken symmetry case or spinodal instabilities 
in the broken symmetry phase. 
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The expression for the renormalized energy density given by ( p.65| - pT67D differs from the 
effective potential in several fundamental aspects: i) it is always real as opposed to the 
effective potential that becomes complex in the spinodal region, ii) it accounts for particle 
production and time dependent phenomena. 

The effective potential is a useless tool to study the dynamics precisely because it misses 
the profuse particle production associated with these dynamical, non-equilibrium and non- 
perturbative processes. 



III. THE UNBROKEN SYMMETRY CASE 



A. Analytic Results 



In this section we turn to the analytic treatment of equations (|2.36|j237|J2.42| ) in the 
unbroken symmetry case. Our approximations will only be valid in the weak coupling 
regime and for times small enough so that the quantum fluctuations, i.e. <?£(t) are not 
large compared to the "tree level" quantities. We will see that this encompasses the times 
in which most of the interesting physics occurs. 

Since S(0) = 0, the back-reaction term g^(r) is expected to be small for small g during 
an interval say < r < T\. This time r 1; to be determined below, determines the relevant 
time scale for preheating and will be called the preheating time. 

During the interval of time in which the back-reaction term gS(r) can be neglected, we 
can solve eq.( |2.36|) in terms of elliptic functions, with the result: 



r](r) = 7]o cn ( T 



k 



^(1 + ril) 



(3.1) 



where cn stands for the Jacobi cosine. Notice that r](r) has period 4u = 4K(k)/^Jl + rj^, 
where K(k) is the complete elliptic integral of first kind. In addition we note that since 



t](t + 2lu) = —7]{t) 



(3.2) 



if we neglect the back-reaction in the mode equations, the 'potential' (— 1 — i] 2 (t)) is periodic 
with period 2a>. Inserting this form for r/(r) in eq.( |2.37D and neglecting (?S(r) yields 



dr 2 



+ q 2 + l + r] 2 - rfi sn 2 rJl + V 2 , k 



Vq {r) = . 



(3.3) 



This is the Lame equation for a particular value of the coefficients that make it solvable in 
terms of Jacobi functions |T7|]. We summarize here the results for the mode functions. The 
derivations are given in the Appendix. 

Since the coefficients of eq.( |3.3| ) are periodic with period 2a; , the mode functions can be 
chosen to be quasi-periodic (Floquet type) with quasi-period 2uj. 



U q (r 



2co) 



e lFM U q (r), 



(3.4) 
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where the Floquet indices F(q) are independent of r. In the allowed zones, F(q) is a 
real number and the functions are bounded with a constant maximum amplitude. In the 
forbidden zones F(q) has a non-zero imaginary part and the amplitude of the solutions either 
grows or decreases exponentially. 

Obviously, the Floquet modes U q [r) cannot obey in general the initial conditions given 
by ( p.20| ) and the proper mode functions with these initial conditions will be obtained as 
linear combinations of the Floquet solutions. We normalize the Floquet solutions as 

17,(0) = 1 • (3.5) 

We choose U q (r) and U q (—r) as an independent set of solutions of the second order differ- 
ential equation ( |3.3[) . It follows from eq. (|3.4j ) that U q (—r) has —F(q) as its Floquet index. 

We can now express the modes f q (r) with the proper boundary conditions (see |2.20| ) as 
the following linear combinations of U q {r) and U q (—r) 



2jn q l 



(3.6) 



where W q is the Wronskian of the two Floquet solutions 

W q = W[U q (r), U q {-r)\ = -217,(0) . (3.7) 



Eq.( |3.3| ) corresponds to a Schrodinger-like equation with a one- zone potential [Tj|. We 
find two allowed bands and two forbidden bands. The allowed bands correspond to 

2 2 

— <g 2 <0 and ^<g 2 <+oo, (3.8) 

and the forbidden bands to 

-oo<g 2 <-l-f and < q 2 < f . (3.9) 

The last forbidden band is for positive q 2 and hence will contribute to the exponential growth 
of the fluctuation function E(r). 

The mode functions can be written explicitly in terms of Jacobi ^-functions for each 
band. We find for the forbidden band, 

TJ I f \ — p-T V~L+vl Z(2K(k) v) ^4(0) ^1(^+2^) m) 

where v is a function of q in the forbidden band < q < defined by 

q=^=cn(2K(k)v,k), 0<v<^. (3.11) 
and Z(u) is the Jacobi zeta function [IS]. It can be expanded in series as follows 

00 ~n 

2 K(k) Z{2K(k) v) = Atc ^ - — sin(2r27ri;) (3.12) 

— 1 1 Q 



n=l 
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where q = e 7rK ( k )/ K ( k > m The Jacobi ^-functions can be expanded in series as follows [20 

oo 

= 2^ (-l) n+1 gC"- 1 ^) 2 S i n ( 2ri - 1)™ , 



ra=l 

00 



i? 4 («|g) = 1 + 2 ^ (-l) n g n cos(2n7ri>) . (3.13) 



71=1 



We explicitly see in eq. (|3.10| ) that U q {r) factorizes into a real exponential with an expo- 
nent linear in r and an antiperiodic function of r with period 2uj. Recall that 

■d x {x + 1) = , $ A (x+l) = +$ 4 (x) . (3.14) 

We see that the solution U q {r) decreases with r. The other independent solution U q (—r) 
grows with r. 

The Floquet indices can be read comparing eq. (|3~^) , ( |3.10j ) and ( |3.14| ), 

F(q) = 2i K(k) Z(2K(k) v) ± tt . (3.15) 
U q {r) turns out to be a real function in the forbidden band. It has real zeroes at 

t = 2uj{n — v) , n eZ . (3.16) 

and complex poles at 

r = 2umi + (2^2 + 1)uj' , ni,n2 eZ . (3-17) 

where u' is the complex period of the Jacobi functions. Notice that the pole positions 
are g-independent, and that U q {r) becomes an antiperiodic function on the borders of this 
forbidden band, q = and q — We find using eq. ( |3.10| ) and ref. fl~9|l , 



U g (r)\ q=0 = cn{r^/l + ^ ,k) (3.18) 
lim [vU q {r)\ = sn(r v /TT^, k) , (3.19) 

respectively. 

The functions U q (r) transform under complex conjugation in the forbidden band as 

[U q (r)}* = U q (r) . (3.20) 
For the allowed band < q < 00, we find for the mode functions 



Uq {r) = e-^ (, ^ o) 1 ' ' ,f (fc) ^ , (3.21) 



where 



g= + (2^)^,^) , (3.22) 
v sn 

0<i;<i , oc > g > (3.23) 
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We see that U q (r) in this allowed band factorizes into a phase proportional to r and a 
complex periodic function with period 2uo. This function U q {r) has no real zeroes in r 
except when q is at the lower border q = . Its poles in r are g-independent and they are 
the same as those in the forbidden band [see eg. ( |3.17| )] . 

The Floquet indices can be read off by comparing eq. (|3.4|) , Q3.14Q and ( |3.21|) 



These indices are real in the allowed band. 

The functions U q {r) transform under complex conjugation in the allowed band as 

[U q (r)}* = U q (-r) . (3.25) 

Obviously these modes will give contributions to the fluctuation S(t) which are always 
bounded in time and at long times will be subdominant with respect to the contributions 
of the modes in the forbidden band that grow exponentially. 

The form of these functions is rather complicated, and it is useful to find convenient 
approximations of them for calculational convenience. 

The expansion of the "^-functions in powers of q = e~ wK ^l K ^ k "> converges quite rapidly 
in our case. Since < k < l/\/2 [see eq. ( |3.1| )], we have 

< q < e- T = 0.0432139 .... (3.26) 

q can be computed with high precision from the series |HJ 

q = A + 2 A 5 + 15 A 9 + 150 A 13 + 1707 A 17 + . . . , 

where (not to be confused with the coupling constant) 

We find from eq. (|3.1|) 



_ 1 (1 + r/Q 2 ) 1 / 4 - (1 + ^/2)^ 

2 (1 + ^) V4 + (1 + ^/2)1/4- ^ 

The quantity A can be computed and is a small number: for < r/ < oo, we find < A < 
0.0432136. . .. Therefore, to very good approximation, with an error smaller than ~ 10~ 7 , 
we may use: 

1 (l + ^o 2 ) 1/4 - (l + ^o 2 /2) 1/4 ,„„ q . 
Q 2 (l+% 2 )V4 + (i + 7/ 2/2)i/4- ^ 

We find in the forbidden band from eq.( |3.10|) and [|T^] 

U ( r ) = e -^T v/l+^O 9 sin(27™) [l+2q (cos 2^-2)+0(g 2 )] 1 ~ 2 Q Sin7r(f + ^j) r + 2 

9 1 — 2q cos(— ) smnv L 

(3.30) 
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where now we can relate v to q in the simpler form 



Vo 

q = —j= COS TTV 

v2 



O 0/ O \ . ^ , q 

1 — Aq sin nv + Aq sin wv (1 + 4 cos 7ri>) + 0(g 



which makes it more convenient to write g(v ) in the integrals, and 



7T 

2~UJ 



1 + r?g 1 - Aq + 12<f + 0(g 4 ) 



(3.31) 



^3.32) 



where < v < \. 

The Floquet indices can now be written in a very compact form amenable for analytical 
estimates 



F(q) = Ain q sin(27ru) 



2q cos2vro + 0(q 2 



+ 71 . 



(3.33) 



In this approximation the zero mode (|3.1|) becomes 



7](t) = 1] COS 



7TT 



2lu 



1 — Aq sin 



,TCT , 



'2a;' 



+ 0(<f) 



(3.34) 



This expression is very illuminating, because we find that a Mathieu equation approxi- 
mation, based on the first term of eq. ( |3.34j ) to the evolution of the mode functions is never 
a good approximation. The reason for this is that the second and higher order terms are 
of the same order as the secular terms in the solution which after resummation lead to the 
identification of the unstable bands. In fact, whereas the Mathieu equation has infinitely 
many forbidden bands, the exact equation has only one forbidden band. Even for small q, 
the Mathieu equation is not a good approximation to the Lame equation p9| . 

From eq. (|3.21|) analogous formulae can be obtained for the allowed band 



U q {r) 



where 



"1? coth 



TV K'(k) 

Wit) 



1 - 2q 



1 - 2q cos(^) 



1 — 2q cos — — 2iv log q 



2q cosh (2 v log q) 



l + 0(,f) , (3.35) 



2 3 / 2 sinhf^f^ 




{l + 2q cosh(2w logg) + 0(q 2 )} . (3.36) 



Here, 



0<v<- 



Note that eq. ( |3.32| ) holds in all bands. 

We can now estimate the size and growth of the quantum fluctuations, at least for 
relatively short times and weak couplings. For small times < r < T\ (to be determined 
consistently later) and small coupling g « 1, we can safely neglect the back-reaction term 
g£(r) in eq.( |2.37| ) and express the modes (p g (r) in terms of the functions U q (r) and U q (—r) 
for this however, we need the Wronskian, which in the forbidden band is found to be given 
by: 
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1 d i #iO>) 

W * = — "T log TTT 



sn 



(3.37) 



In terms of the variable q 2 this becomes, after using eq. ( |3. 1 1| ) 



-2g 



\ 



+ 1 + g 2 



2 



(3.38) 



This Wronskian is regular and non-zero except at the four borders of the bands. 
We find from eq.fl3.6p that |<^ 9 (t)| 2 is given by 



1 f 4Q 2 

<P,(r) | 2 = ^ | [U q {r) + U q (-r)] 2 + [17 s (t) - U q {-r)\ 



^3.39) 



where we took into account eqs. ( |3.20| and ( |3.25| ). Notice that both terms in the rhs of 



eq.( 3.39 ) are real and positive for real q. For very weak coupling and after renormalization, 
the contribution to g"E(r) from the stable bands will always be perturbatively small, while the 
contribution from the modes in the unstable band will grow in time exponentially, eventually 
yielding a non-perturbatively large contribution. Thus these are the only important modes 
for the fluctuations and the back-reaction. An estimate of the preheating time scale can be 
obtained by looking for the time when ^S(r) is of the same order of the classical contributions 
to the equations of motion. In order to obtain an estimate for the latter, we consider the 
average over a period of the classical zero mode: 



l+<V 2 (r) >=(l+v 2 o) 



2E(k) 
K{k) 



- 1 



(3.40) 



which yields for small and large initial amplitudes the following results 



< V 2 (r) >"°=° | 
<r/ 2 ( r ) >^°° 0.4569 



(3.41) 
(3.42) 



Therefore the average over a period of r] 2 (r) is to a very good approximation i] 2 /2 for all 
initial amplitudes. This result provides an estimate for the preheating time scale 7~i; this 
occurs when gE(ri) ~ (1 + r] 2 /2). Furthermore, at long times (but before gE (1 + rj 2 /2)) 
we need only keep the exponentially growing modes and <?E(r) can be approximated by 



4 Jo 



9 fV2 2 , - 
q aq — 



1 



i L 



4^1 
W 2 qi 



U q {~r) 



(3.43) 



Moreover, choosing r such that the oscillatory factors in U q (—r) attain the value 1 (the 
envelope), and using eq. fl3.10|) we finally obtain: 



J est—env 



1 f— I 

;r) = 4Jo 2q2dq n a 



w 2 qi 



,2r yjl+nl Z(2K(k)v,k) 



(3.44) 
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where v depends on the integration variable through eq. (|3.11|) . 

The Jacobi Z function can be accurately represented using eq. (|3.12| ) 

Z{2K(k) v,k) =4q sin 2ttv [1 - 2 q (2 - cos 2ttv)} + 0(q 3 ) . 



(3.45) 



where we recall that q < 0.0433. 

The integral Q3.44| ) will be dominated by the point q that maximizes the coefficient of r 
in the exponent. This happens at q = g 1; v — v ii where 



?i = \vo (1 - q) + 0{q 2 ) 
Z(2K(k) v u k) = Aq (1 - 4g) + 0(<f) 

We can compute the integral ( 3.44|) by saddle point approximation to find: 



(3.46) 
(3.47) 



q\ 



402 
^ W 2 



J est—env \ 



2 a 



dq e -^ T{q - qi)2 * v 71 ^ d- 2 « [1 + 0(g)] 

e 8r v /i+^g 3(1-4$) 



4Q 2 
1 + 1% 



W 2 



l + 0( 



64(1 + % 2 )V4 

We can relate g to ?7o using eq.( 3.29|) , and we have used the small q expansion 

Z"(2K(k) v u k) = -lQq (1 - Aq) + 0(f) 

9q) + 0(q 2 ). 



dv 



-5P 



(3.48) 

(3.49) 
(3.50) 



In summary, during the preheating time where parametric resonance is important, 
^est-envij) can be represented to a very good approximation by the formula 



->est—env \ 



(3.51) 



where B and N are functions of 770 given by 



B = 8^1+ V 2 q(l-4q) + O(q 3 ) 
64 (l+r/o 2 ) 1 ^ ^n qi 



N 



71 



1/2 



40 2 
1 + 1% 



w 2 



(4 + 3^)^4 + 5^ 



% 3 (1 + % 2 ) 3/4 



(3.52) 



and eq. ( |3.29| ) gives q as a function of t]q . This is one of the main results of this work. 
We display in Table I below some relevant values of q, B and N as functions of t/q. 
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We notice that the limiting values of B and N for 770 
mation even for 770 ~ 1. Namely, 



/ »n3 ,,-Boo Vo t 



with the asymptotic values given by 

Bnc 



8e- n (l-4e~ n )[l + 0(^ 
12 - -*/2r 



-=V5 e-*/ 2 [l + 0(% 2 )] = 3.147 ... [1 + O(r^)]. 

'7T 



00 yield a very good approxi- 

(3.53) 

(3.54) 
(3.55) 



0.285953... [1 + O(r)o 2 )} 

2\ 



These rather simple expressions ( p. 51 - 3*35 ) allow us to perform analytic estimates with 
great accuracy and constitute one of our main analytic results. The accuracy of this result 
will be discussed below in connection with the full numerical analysis including back-reaction. 

Using this estimate for the back-reaction term, we can now estimate the value of the 
preheating time scale T\ at which the back-reaction becomes comparable to the classical 
terms in the differential equations. Such a time is defined by <?£(ti) ~ (1 + tJq/2). From the 
results presented above, we find 



1 N(l + 4/2) 



(3.56) 



The time interval from r = to r ~ n is when most of the particle production takes 
place. After r ~ n the quantum fluctuation become large enough to shut-off the growth of 
the modes and particle production essentially stops. We will compare these results to our 
numerical analysis below. 

We can now use our analytic results to study the different contributions to the energy and 
pressure coming from the zero mode and the quantum fluctuations and begin by analyzing 
the contribution to the energy e and pressure p from the zero mode rj(r). 

The dimensionless energy and pressure, (normalized by the factor 2M^/ A#) are given by 
the following expressions, 



Po(r) 



V 2 + V{r) 2 + -v(r) A 
V 2 ~ V(r) 2 ~ ^( r ) 4 



When the back-reaction term <?£(t) can be neglected, we can use eq. 
approximation to T](t). In this approximation 



(3.57) 



T[) as a good 



e 



1 + 2*8 



Po(r) + 



e = r,l (l + r,t) sn 2 dn 2 (ry/l + ^, k^j . 



(3.58) 



The zero mode energy is conserved and the pressure oscillates between plus and minus eo 
with period 2u. 
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Averaging poij) over one period yields 



2*c 



2lu Jo 

Inserting eq. (|3.58|) into eq. (|3.59| ) yields |24| 



dr po(r) 



(3.59) 



1 2 



i - 



(3.60) 



where A; is given by eq. (|3.1|) . 

< po > vanishes for small rjo faster than eo, 



(3.61) 



so that the zero mode contribution to the equation of state is that of dust for small r/ . For 
large rj we find from eq. ( |3.60 ), 



rio^oo 1 4 2 

< Po > = + % 



2 £(lA/2) 



+ 0(1) 



(3.62) 



w here ± - 2 g(1/v ^ 

wnele 2 3 Jf(l/V3) 

r/ -> oo: 



0.01437. 



2 3 K(l/V2)_ 
The equation of state approaches that of radiation for 



J70->oo 

< po > = e 



0.6092... 



(3.63) 



Thus we see that for small amplitudes the zero mode stress-energy, averaged over an 
oscillation period, behaves as dust while for large amplitudes, the behavior is that of a 
radiation fluid. The ratio < p > /e for zero mode vs.e is shown in figure 1. 

The contribution from the k ^ modes originates in the quantum fluctuations during 
the the stage of parametric amplification. 

Since we have fluid behaviour, we can define an effective (time-dependent) polytropic 
index 7(7") as 



7(r) 5 ^ + l. 



(3.64) 



where renormalized quantities are understood throughout. Within a cosmo logical setting 
whenever 7(7") reaches a constant value such equation of state implies a scale factor R(r) = 

In the case being studied here, that of Minkowski space, e is time-independent and hence 
equal to the initial energy density (divided by N and restoring pre-factors) which after a 
suitable choice of the constant C is given by: 



2IM 



R 



4 '1 2 

2 Vo 



1 + 2* 



(3.65) 
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As argued before, for weak coupling the important contribution to the quantum fluc- 
tuations come from the modes in unstable bands, since these grow exponentially in time 
and give rise to a non-perturbatively large contribution. Thus we concentrate only on these 
modes in calculating the pressure. 

The contribution of the forbidden band to the renormalized p(r) + e can be written as 



\p(r)+s[ 



2\Mr 



unst 



A 



9 



R 



V0/V2 



q 2 dq 



( Pq\ T ) 



1 2 1 2 

1 + 3 q 



(3.66) 



After renormalization, the terms that we have neglected in this approximation are pertur- 
batively small (of order g) whereas the terms inside the bracket eventually become of order 
1 (comparable to the tree -level contribution). We now only keep the exponentially growing 
pieces in the mode functions f q (T) an d 0<j(t) since these will dominate the contribution to 
the pressure. This is simplified considerably by writing to leading order in q 



1 



?7o cot 



Averaging over a period of oscillation yields 



TC(V 



2a/ 



+ 0{q) 



2\M, 



R 



\unst 



9 
4 



V2 q dq 



(4vr) 2 2 n q sm 2 nv 



1 H 1 

W 2 



d 2t ^Jl+r,l Z{2K(k)v,k) 



(3.67) 

(3.68) 
(3.69) 



This integral is similar to the one in eq.( |3.44| ) and we find that they are proportional in 
the saddle point approximation. In fact, 



\p{r) + e] unst 



2\M 



R\ 



gz 



est— env 



T 



1 13 2 
1+ 12"» 



(3.70) 



where T, est _ env (r) is given by eq. (ft.51|) . 
The effective polytropic index j(t) is: 



7(r) = gT, 



est— env 



T 



12 + 13^ 



(3.71) 



3*0 + 2) ' 

When 5 f S es t-eni)(' r i) ~ 1 + t/q/2, i.e. at the end of the preheating phase, j(t) is given by 

(3.72) 



leff OC 



12 + 13r7g 



6r/ 2 



We note here that for very large 770 the effective polytropic index is 7 e // — 13/6 ~ 0{1). It 
is clear then that the physics can be interpreted in terms of two fluids, one the contribution 
from the zero mode and the other from the fluctuations, each with an equation of state that 
is neither that of dust nor of radiation, but described in terms of an effective polytropic 
index. 

We can now use our approximations to obtain an estimate for the number of particles 
produced during the preheating stage. In terms of dimensionless quantities, the particle 
number, defined with respect to the initial Fock vacuum state is given by eq.( [2.53| ). 
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This particle number will only obtain a significant contribution from the unstable modes 
in the forbidden band where to leading order in q we can approximate v ? g( r ) an d ^(t) by 
its exponentially growing pieces [see eq. Q3.39Q 1, as follows: 



\<Pq\T) 



4 a 



4 ft 2 
1 H 1 

w 2 



U q {~r) 



1 + r]l) cotg 2 



7T V — 



2u' 



Mr)\ 2 + 0(q) . 



(3.73) 
(3.74) 



The total number of produced particles N(t) per volume |M#| 3 is given by: 



N{t) 
\M R \- 



d 3 q 



N q {r) . 



(3.75) 



The asymptotic behaviour (|2.26|) ensures that this integral converges. 
Following the same steps as in eq. (|3.44 ) and ( |3.68| ), we find 



N{t) 



wrist 



8tt 2 



J est—env 



qi 



1 4+| % 2 



A/ 



R V 4 + 5 ^o 



est— env 



(3.76) 



where we used eq. (|3.46| ) and S est _ ent) (r) is given by the simple formula ( |3.51| ). Notice that 
by the end of the preheating stage, when (?S(r) m 1 + i]q/2 the total number of particles 
produced is non-perturbatively large, both in the amplitude as well as in the coupling 



tot 



1 (4+1^(1 + ^/2) 
X R y/4 + 5 7? 2 



(3.77) 



The total number of adiabatic particles can also be computed in a similar manner with 
a very similar result insofar as the non-perturbative form in terms of coupling and initial 
amplitude. 



B. Numerical Results 

We now evolve our equations for the zero and non-zero modes numerically, including the 
effects of back-reaction. We will see that up to the preheating time, our analytic results 
agree extremely well with the full numerical evolution. 

The procedure used was to solve equations (|2.36j , |2.37|) with the initial conditions (|2.38| , 



2.39 , 2.40|) and (|2.42|) using a fourth order Runge-Kutta algorithm for the differential equation 



and an 11-point Newton-Cotes integrator to compute the fluctuation integrals. We tested 
the cutoff sensitivity by running our code for cutoffs A/|Mr| = 100, 70, 50, 20 and for very 
small couplings (which is the case of interest. We found no appreciable cutoff dependence. 
The typical numerical error both in the differential equations and the integrals are less than 
one part in 10 9 . 

Figure 2. a shows i](t) vs.r for r] = 4.0 ,g = 10~ 12 . For this weak coupling, the effect 
of back-reaction is negligible for a long time, allowing several undamped oscillations of the 
zero mode. Figure 2.b shows (?£(t) vs. r. It can be seen that the back-reaction becomes 
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important when ^S(t) ~ 1 + i] 2 /2 as the evolution of 77 (r) begins to damp out. This 
happens for r ~ 25 in excellent agreement with the analytic prediction given by eq.( p.56 ) 
7~i = 26.2 . . ., the difference between the analytic estimate for E(r) given by eq.( |3.53| ) and the 
numerical result is less than 5% in the range < r < 30. Figure 2.c shows gj\f{r) vs. r and 
we see that the analytic expression ( |3.77| ) gives an approximate estimation XnAftot ~ 74.6 . . . 
for the final number of produced particles. 

Figures 2.d-2.f, show gN q {r) for r = 40, 120, 200, we see that the prediction of the width 
of the unstable band < q < r/o/v^2 is excellent and is valid even for very long times beyond 
the regime of validity of the small time, weak coupling approximation. However, we see that 
the peak becomes higher, narrower and moves towards q ~ 0.5 as time evolves beyond t%. 
This feature persists in all numerical studies of the unbroken phase that we have carried 
out, this changes in the peak width, height and position are clearly a result of back-reaction 
effects. We have searched for unstable bands for < q < 20 and we only found one band 
precisely in the region predicted by the analytic estimate. All throughout the evolution there 
is only one unstable band. The band develops some structure with the height, position and 
width of the peak varying at long times but no other unstable bands develop and the width 
of the band remains constant. For values of q outside the unstable band we find typically 
gN q < 10~ 13 at all times. This is a remarkable and unexpected feature. 

Obviously this is very different from the band structure of a Mathieu equation. The 
Mathieu equation gives rise to an infinite number of narrowing bands, so that quanti- 
tative estimates of particle production, etc. using the Mathieu equation approximation 
would be gross misrepresentations of the actual dynamics, with discrepancies that are non- 
perturbatively large when the back- reaction becomes important |29j . Since particle produc- 
tion essentially happens in the forbidden bands, the quantitative predictions obtained from 
a single forbbiden band and an infinite number, as predicted by WKB or Mathieu equation 
analysis, will yield different physics. 

We have carried the numerical evolution including only the wave-vectors in the unstable 
region and we find that this region of g-wavevectors is the most relevant for the numer- 
ics. Even using a cutoff as low as q c = 4 in this case gives results that are numerically 
indistinguishable from those obtained with much larger cutoffs q c = 70 — 100. The occu- 
pation number of modes outside the unstable bands very quickly becomes negligible small 
and for q m 4 it is already of the same order of magnitude as the numerical error < 10~ 10 . 
Clearly this is a feature of the weak coupling case under consideration. Keeping only the 
contribution of the modes in the unstable band, the energy and pressure can be written as 
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(3.78) 
(3.79) 

(3.80) 



where we have made explicit that we have neglected terms of order g in eqs.( |3.78| - |3.79| ). The 
terms multiplied by g in eqs. ( [3.78| , |3.79|) become of order 1 during the preheating stage. For 
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the parameters used in figures 2, we have checked numerically that the energy Q3.78| ) is con- 
served to order g within our numerical error. Figure 2.g shows the pressure 2|M/j| 4 p{t)/\r 
versus r. Initially p(0) = —e (vacuum dominated) but at the end of preheating the equation 
of state becomes almost that of radiation p^ = e/3. 

For very small coupling (g ~ 1CT 12 ), the backreaction shuts-off suddenly the particle 
production at the end of the preheating (see fig. 2c). Later on, (r larger than 100 for 
g ~ 10~ 12 ) the time evolution is periodic in a very good approximation. That is, this non- 
linear system exhibits a limiting cycle behaviour. The modulus of the fc-modes do not grow 
in time and no particle production takes place. This tells us that no forbidden bands are 
present for q 2 > in the late time regime. 

We have numerically studied several different values of t]o ,g finding the same qualitative 
behavior for the evolution of the zero mode, particle production and pressure. In all cases 
we have found remarkable agreement (at most 5% difference) with the analytical predictions 
in the time regime for which < gS(r) < 1. The asymptotic value of the pressure, however, 
only becomes consistent with a radiation dominated case for large initial amplitudes. For 
smaller amplitudes rjo = 1 we find that asymptotically the polytropic index is smaller than 
4/3. This asymptotic behavior is beyond the regime of validity of the approximations in 
the analytic treatment and must be studied numerically. This polytropic index depends 
crucially on the band structure because most of the contribution comes from the unstable 
modes. 



IV. THE BROKEN SYMMETRY CASE 



A. Analytic Results 



As in the unbroken case, for g << 1 we can neglect gS(r) in eq. (|2.43 ) until a time r 2 
at which point the fluctuations have grown to be comparable to the 'tree level' terms. 
We then find for < r/o < 1, 

Vo 



7](T) 



dn(r\/l- f 



l-V 2 o 



(4.1) 



Notice that r\ (r) has period 2uj 
For 1 < ?7o < \f2 we find 



2K{k) 



1 2 



. The elliptic modulus k is given by eq.(|4.1|). 



T)(t) = r] dn (r^o/V^fe 



2(1 



(4.2) 

This solution follows by shifting eq.( [4.1| ) by a half-period and changing rj^ — > 2 — rj 2 . It has 
a period 2u = ^^-K(k). For i] — > 1, 2u — > n^/2 and the oscillation amplitude vanishes, 
since rj = 1 is a minimum of the classical potential. 
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For 7/0 > Vz we obtain 



t/(t) = 77 cn (yT/g - lr, A; 
?7o 



fc 



2(r/o 2 



(4.3) 



This solution has 4w = 4 ^j_ as period. 

The solutions for r?o < v^2 and r?o > y/2 are qualitatively different since in the second case 
t/(t) oscillates over the two minima 7/ = ±1. In the limiting case 77 = \f2 these solutions 
degenerate into the instanton solution 



r/(r) 



y/2 
coshr 



(4.4) 



and the period becomes infinite. 

Inserting this form for 77 (t) in eq.( 2.44j) and neglecting #S(t) yields for < tjq < 1, 



dr 2 



+ q 2 -l + 



dn 2 (r\/l-f,fc 



ip q {r) = 0. 



(4.5) 



This is again a Lame equation for a one- zone potential and can also be solved in closed form 
in terms of Jacobi functions. We summarize here the results for the mode functions, with 
the derivations again given in the Appendix. 

As for unbroken symmetry case, there are two allowed bands and two forbidden bands 
The allowed bands for < r? < 1 correspond to 



2 

< <7 2 < — and 1 
- H ~ 2 



I < q 2 < +00 



and the forbidden bands to 



2 2 

00 < q 2 < and ^ < q 2 < 1 - ^ 
-h - 2 -h - 2 



(4.6) 



(4.7) 



The last forbidden band exists for positive q 2 and hence contributes to the growth of £(t). 

The Floquet solutions obey eqs. ( |3.4|) and (|3~5| ) and the modes <p q (r) can be expressed 
in terms of U q (r) and U q (—r) following eq.([ 



It is useful to write the solution U q (r) in terms of Jacobi -^-functions. For the forbidden 

2 2 

band ^ < q 2 < 1 — after some calculation (see the Appendix), 
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where < v < h is related with q through 
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and k is a function of i] as denned by eq.( |4.1| ). 

We see explicitly here that U q {r) factorizes into a real exponential with an exponent 
linear in r and an antiperiodic function of r with period 2uo. 

The Floquet indices for this forbidden band are given by 

F(q) = 2iK(k) Z(2K(k)v) ± tt . (4.10) 

2 

For the allowed band 1 — ^ < q 2 < +oo , we find for the modes, 

,l > MB Mi;) • ( ' 

where q and a are related by: 



i _ 2a 

2 



sn ( cn\/l ^J, 



with > a > 0. The Floquet indices for this first allowed band are given by 



2 

Analogous expressions hold in the other allowed band, < q 2 < ^ 



-jr. tf 3 (0) ^ 4 (^) 



^(g) Mis) 



(4.13) 



are given by 



Here, q = ^= sn (ay 1 — A;' j and ^ > a > 0, and the Floquet indices for this band 



$2 / \ 



F M = < sUJ- (414) 



For ?7o ~ 1 the situation is very similar to the unbroken symmetry case; the zero mode 
oscillates quasi-periodically around the minimum of the tree level potential. There are effects 
from the curvature of the potential, but the dynamics can be analyzed in the same manner 
as in the unbroken case, with similar conclusions and will not be repeated here. 

The case 1] << 1 is especially interesting || for broken symmetry because of new 
and interesting phenomena |§ that has been recently associated with symmetry restoration 

In this limit, the elliptic modulus k [see eq. ( |4.1|) 1 approaches unity and the (real) period 
2ui grows as 



2 oo ~ 2 K(k) + O( V 2 ) ~ 2 In f ^ j + 0( V 2 Q ) (4.15) 
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In this limit, both the potential in eq. (|4.5|) and the mode functions ( |4.8|) can be approximated 



by hyperbolic functions [19]: 



where 
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< u < A 



Using the imaginary Jacobi transformation [19| 
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we can derive expressions for the mode functions U q {r) valid for small 770: 



U q (r) 



-t tanh«(l+^h 1 — % COshM COsh(« 
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(4.22) 



Here u is related with q through eq. Q4.18|) . We see that the function U, 
almost as e T for q near the lower border of the forbidden band u 
can be interpreted as the joint effect of the non-periodic exponential factor in eq 
the growth of the periodic ^-functions. Since the real period is here of the order A, the 
two effects cannot be separated. The unstable growth for r < uj also reflects the spinodal 



,., r) grows with r 
A. This fast growth 
and 



instabilities associated with phase separation 11]. 

In this case, there is a range of parameters for which the quantum fluctuations grow 
to become comparable to the tree level contribution within just one or very few periods. 
The expression ( (4.22j ) determines that S(r) m e 2r from the contributions of modes near the 
lower edge of the band. The condition for the quantum fluctuations to become of order 
1 within just one period of the elliptic function is ge iuJ ~ 1 which leads to the conclusion 
that for 77(0) < g 1 / 4 the quantum fluctuations grow very large before the zero mode can 
actually execute a single oscillation. In such a situation an analysis in terms of Floquet 
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(quasi-periodic) solutions is not correct because the back reaction prevents the zero mode 
from oscillating enough times for periodicity to be a reasonable approximation. 

We now analyze the behavior of the pressure for the zero mode to compare to the previous 
case. In the approximation where eqs.( f4.1|) , ( }4.2|) and (|4.3| ) hold and adjusting the constant 
C in the definition of the energy, we have 



e 



P0{T) 



-e + v(t) 



(4.23) 



Inserting eqs.(|4~l"l), ( ^2|) and (fO| ) in eq.( |4.23|) yields 
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% > V2 : 
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Notice that the functional form of the elliptic modulus k as a function of rjo is different in 
each interval [see eqs.Q, flOD and 

Let us now average the pressure over a period as in eq. (|3.59|) . We find 
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The dimensionless energy eo tends to 1/4 both as 770 — > and 7/0 — > V% in both cases we 
find using eq. ( |4.25| ) 
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This is result is recognized as vacuum behaviour in this limit. 
For 770 — > 1, eq. fl4.25| ) yields, 



< Po > vo^i 



O(vo 



V 



(4.27) 



That is a dust type behaviour, which is consistent with the small amplitude limit of the 
unbroken symmetry case studied before. 

Finally, for r] —>■ 00, when the zero mode is released from high up the potential hill, we 
find that the pressure approaches radiation behaviour (from above) 
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Figure 3 shows < po > /eq vs. Eq. 

As mentioned before, we expect that for t)q « 1 the conclusion will be modified dra- 
matically by the quantum corrections. 



B. Numerical Results 

In the region t]q m 1 the analytic estimates are a good approximation for large times and 
weak couplings. We have studied numerically many different cases with rjo > 0.5 and weak 
coupling and confirmed the validity of the analytic estimates. These cases are qualitatively 
similar to the unbroken symmetry case with almost undamped oscillations for a long time 
compatible with the weak coupling approximation and when gT^{r) grows by parametric 
amplification to be of order one with a consequently large number of produced particles and 
the evolution of the zero mode damps out. 

However as argued above, for tjq << 1 the analytic approximation will not be very 
reliable because the quantum fluctuations grow on a time scale of a period or so (depending 
on the coupling) and the back-reaction term cannot be neglected. Thus this region needs to 
be studied numerically. 

We numerically solved equations (|2.43|) - (|2.44j) with the initial conditions ( 2.42|) , Q2.45|) , 



( j2.46| ) and ( |2.47| ). The numerical routines are the same as in the unbroken symmetry case. 
Again we tested cutoffs A/|M R | = 100, 70, 50, 20 and for very small couplings (which is the 
case of interest, g = 10~ 6 ---g = 10~ 12 ) we found no appreciable cutoff dependence, with 
results that are numerically indistinguishable even for cutoffs as small as q c ~ 2 . The typical 
numerical error both in the differential equations and the integrals are the same as in the 
unbroken case, less than one part in 10 9 . 

We begin the numerical study by considering first the case of very small coupling and 
rjo << 1; later we will deal with the case of larger couplings and initial values of the zero 
mode. Figure 4. a shows i](r) vs. r for r/ = 10~ 5 , g = 10~ 12 . In this case we see that 
within one period of the classical evolution of the zero mode, ^S(t) becomes of order one, 
the quantum fluctuations become non-perturbatively large and the approximation valid for 
early times and weak couplings breaks down. Fig. 4.b shows ^S(r) and fig. 4.c shows gj\f{r) 
vs. r for these parameters. We find that only the wave vectors in the region < q < 1 are 
important i.e. there is only one unstable band whose width remains constant in time. This is 
seen in figs. 4.d-f, which show the particle number (defined with respect to the initial state) 
as a function of wave vector for different times, gN q {r = 30) ,gN q (r = 90) ,gN q (r = 150) 
respectively. Although the analytic approximation breaks down, the prediction eq. ( j4.7| ) for 
the band width agrees remarkably well with the numerical result. As in the unbroken case, 
the band develops structure but its width is constant throughout the evolution. As can 
be seen in these figures the peak of the distribution becomes higher, narrower and moves 
towards smaller values of q. The concentration of particles at very low momentum is a 
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consequence of the excitations being effectively massless in the broken symmetry case. The 
features are very distinct from the unbroken symmetry case, in which the peak approaches 
q pa 0.5. 

We found in all cases that the asymptotic behavior corresponds to 



M 2 (r) = -l+ V 2 (T)+gE(r: 



lirm 







(4.29) 



This is a consistent asymptotic solution that describes massless "pions" and broken symme- 
try in the case 77(00) 7^ 0. 

For times r pa 100 — 150 the value of the zero mode is somewhat larger than the initial 
value: 77 (r = 150) pa 2 x 10~ 5 . This result, when combined with the result that the average 
of the effective mass approaches zero is clearly an indication that the symmetry is broken. 
We found numerically that the final value of the zero mode depends on the initial value and 
the coupling and we will provide numerical evidence for this behavior below. 

Figure 4. a, presents a puzzle. Since the zero mode begins very close to the origin with 
zero derivative and ends up very close to the origin with zero derivative, the classical energy 
of the zero mode is conserved. At the same time, however, the dynamical evolution results in 
copious particle production as can be seen from fig. 4.c. We have shown in a previous section 
that the total energy is conserved and this was numerically checked within the numerical 
error. Thus the puzzle: how is it possible to conserve the total energy, conserve the classical 
zero mode energy and at the same time create 0(l/g) particles? The answer is that there is 
a new term in the total energy that acts as a 'zero point energy' that diminishes during the 
evolution and thus maintains total energy conservation with particle production. The most 
important contribution to the energy arises from the zero mode and the unstable modes 
< q < q u . The energy and pressure are given by (adjusting the constant C such that the 
energy coincides with the classical value): 
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(4.30) 

(4.31) 
(4.32) 
(4.33) 

(4.34) 
(4.35) 



where A4(r) 2 is the effective squared mass of the N — 1 "pions" and again O(g) stand for 
perturbatively small terms of order g. The terms displayed in (|4.31[ - f4.34|) are all of 0(1) 
during the preheating stage. 

We find that whereas En grows with time, the term eq becomes negative and decreases. 
In all the cases that we studied, the effective mass j\4(r) approaches zero asymptotically; 
this is seen in figure 4.g for the same values of the parameters as in figs. 4.a-c. This behavior 
and an asymptotic value rj^ 7^ is consistent with broken symmetry and massless pions by 
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Goldstone's theorem. The term ec in eg. (|4.33|) can be identified with the 'zero' of energy. 
It contributes to the equation of state as a vacuum contribution, that is pc = —£c an d 
becomes negative in the broken symmetry state. It is this term that compensates for the 
contribution to the energy from particle production. 

This situation is generic for the cases of interest for which tjq << 1, such is the case 
for the slow roll scenario in inflationary cosmology. Figs. 4.h-k show e c i(r) vs. r, £at(t) 
vs. r, ec(t) vs. r and 2jM^p -P( r ) vs r ^ or ^ ne same values of parameters as fig. 4. a. 
The pressure has a remarkable behavior. It begins with p = — e corresponding to vacuum 
domination and ends asymptotically with a radiation-like equation of state p = e/3. A simple 
explanation for radiation-like behavior would be that the equation of state is dominated by 
the quantum fluctuations which as argued above correspond to massless pions and therefore 
ultrarelativistic. It must me noticed that we obtain a radiation-like equation of state in spite 
of the fact that the the distribution is out of equilibrium and far from thermal as can be 
seen from figs. 4.d-f. 

An important question to address at this point is: why does the zero mode reach an 
asymptotic value different from the minimum of the effective potential? The answer to 
this question is that the effective potential is an irrelevant quantity to study the dynamics 
PJTT[] . Once there is profuse particle production, a feature completely missed by the effective 
potential, the zero mode evolves in a non-equilibrium bath of these excitations. Through the 
time evolution, more of these particles are produced and the zero mode evolves in a highly 
excited, non-equilibrium state. Furthermore we have seen in detail that this mechanism of 
particle production modifies dramatically the zero point origin of energy and therefore the 
minimum of the effective action, which is now the appropriate concept to use. The final 
value reached by the zero mode in the evolution will be determined by all of these non- 
perturbative processes, and only a full numerical study will capture the relevant aspects. 
As we have argued above, approximations based on Mathieu-type equations or the WKB 
approximation are bound to miss important details and will lead to a completely different 
evolution. 



V. SYMMETRY RESTORATION AT PREHEATING? 

The numerical result depicted by figures 4.a-c, which are very similar to results obtained 
previously ||, has motivated the suggestion that the growth of quantum fluctuations is so 
strong that the non-equilibrium fluctuations restore the symmetry @JT^]. The argument is 
that the non-equilibrium fluctuations given by the term X(ijj 2 (t)) in eq. [2.20| for the mode 



functions grow exponentially and eventually this term overcomes the term — |m 2 | leading to 
an effective potential with a positive mass squared for the zero mode (see eq. |2.19|) . 



Although this is a very interesting suggestion, it is not borne out by our numerical 
investigation. The signal for broken or restored symmetry is the final value of the zero 
mode when the system reaches an equilibrium situation. Any argument about symmetry 
restoration based solely on the dynamics of the fluctuation term \(ip 2 (t)) is incomplete if it 
does not address the dynamics of the zero mode. In particular for the case of figures 4.a-c, 
the initial value of the zero mode rjo ^ and the final value is very close to the initial value 
but still different from zero. 



33 



At the same time, the asymptotic effective mass of the 'pions' is on average zero. Clearly 
this is a signal for symmetry breaking. Because the initial and final values of the order 
parameter are so small on the scale depicted in the figures, one is tempted to conclude that 
the symmetry originally broken by a very small value of the order parameter is restored 
asymptotically by the growth of non-equilibrium fluctuations. To settle this issue we show a 
different set of parameters in figures 5.a,b that unambiguously show that the final value of 
the order parameter rj^ ^ 0, while the effective mass of the pions M.(t) — > 0. Here r/o = 0.01 
and g = 10~ 5 , and asymptotically we find r](r = 150) ~ 0.06, the average of the effective 
mass squared J\4 2 (t) =0 and the symmetry is broken, despite the fact that the fluctuations 
have grown exponentially and a number of particles 0(1/ g) has been produced. 

The reason that the symmetry is not restored is that when the effective mass becomes 
positive, the instabilities shut off and the quantum fluctuations become small. When this 
happens gH is no longer of order one and the instabilities appear again, producing the oscil- 
latory behavior that is seen in the figures for <?£(t) at long times, such that the contributions 
of the oscillatory terms average to zero. It is rather straightforward to see that there is a 
self-consistent solution of the equations of motion for the zero mode and the fluctuations 
with constant 7]^ and A4 2 (oo) = 0. Eq. ( |2.43| ) takes the asymptotic form 



. 



In addition, eq. (|2.44j) yields when M(oo) 2 = 0, 

where A q and B q depend on the initial conditions and g. We get from eqs. Q2.42l ) and ( f2.53| ) 



V l = l-4g «^ Nq (oo) dq-gS( Vo ), (5.1) 
Jo q 2 + \V g 



where S( Vo ) = f (1 - i] 2 ) [log ^ - f 1 + |(1 + V 2 ) ArgTh - arctg 
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2 1 

We see that the value of r]^ depends on the initial conditions. Whereas the last term 
in eq. ( |5.1| ) is perturbatively small, the contribution from the produced particles is non- 
perturbatively large, as N q (oo) ~ 1/g for the unstable wavevectors. Thus the asymptotic 
value of the zero mode is drastically modified from the tree level v.e.v (in terms of renor- 
malized parameters) because of the profuse particle production due to the non-equilibrium 
growth of fluctuations. 

Another way to argue that the symmetry is indeed broken in the final state is to realize 
that the distribution of "pions" at late times will be different than the distribution of the 
quanta generated by the fluctuations in the a field, if for no other reason than that the pions 
are asymptotically massless while the a quanta are massive, as long as rj^ is different from 
zero. If the symmetry were restored during preheating, these distributions would have to be 
identical. 

The conclusion is that the final value of the zero mode depends strongly on the ini- 
tial conditions and couplings, and though symmetry restoration can take place for other 
situations. 

A consistent study of the evolution of the zero mode and quantum fluctuations determines 
what happens in each case f29fl. Thus we emphasize that the ultimate test for symmetry 
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restoration is the asymptotic in time evolution of the zero mode, which is the order parameter 
for symmetry breaking. 



VI. THE REHEATING TEMPERATURE 



The arena in which these results become important is that of inflationary cosmology. 
In particular, the process of preheating is of vital importance in understanding how the 
big-bang cosmology is regained at the end of inflation, i.e. the reheating mechanism. 

While our analysis has been entirely a Minkowski space one, we can make some comments 
concerning the reheating temperature. However, a more detailed analysis incorporating the 
expanding universe must eventually be done along the lines suggested in this work, to get 
more accurate results. 

Since the particles created during the preheating stage are far from equilibrium, thermal- 
ization and equilibration will be achieved via collisional relaxation. In the approximation 
that we are studying, however, collisions are absent and the corresponding contributions are 
of 0(1/N) 0. The difficulty with the next order calculation and incorporation of scattering 
terms is that these are non-local in time and very difficult to implement numerically. 

However we can obtain an estimate for the reheating temperature under some reasonable 
assumptions: in the cosmological scenario, if 'the equilibration time is shorter than the inverse 
of the expansion rate H, then there will not be appreciable redshifting of the temperature 
because of the expansion and we can use our Minkowski space results. 

The second assumption is that the time scales between particle production and ther- 
malization and equilibration are well separated. Within the large N approximation this is 
clearly correct because at large N, scattering processes are suppressed by 1/N. If these 
two time scales are widely separated then we can provide a fairly reliable estimate of the 
reheating temperature as follows. 

Equilibration occurs via the redistribution of energy and momentum via elastic collisional 
processes. Assuming that thermalization occurs on time scales larger than that of particle 
production and parametric amplification, then we can assume energy conservation in the 
scattering processes. Although a reliable and quantitative estimate of the reheating temper- 
ature can only be obtained after a detailed study of the collisional processes which depend 
on the interactions, we can provide estimates in two important cases. If the scattering pro- 
cesses do not change chemical equilibrium, that is conserve particle number, the energy per 
particle is conserved. Since the energy (density) stored in the non-equilibrium bath and the 
total number of particles per unit volume shown in the previous sections: 



with proportionality constants of order one, we can estimate the reheating temperature to 
be: 



N 



£ 



Mr 
\r 

Mr 
Xr 



(6.1) 



(6.2) 




(6.3) 
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Here \Mr\ is the inflaton mass. This is consistent with previous results 0. 

This result seems puzzling, since naively one would expect e ~ T^; iV rs Tjj but the 
powers of A do not match. This puzzle arises from intuition based on an ultrarelativistic 
free particle gas. However the "medium" is highly excited with a large density of particles 
and the "in medium" properties of the equilibrated particles may drastically modify this nave 
result as is known to happen in most theories at high temperature, where the medium effects 
are strong and perturbation theory breaks down requiring hard thermal loop resummation. 

In the case in which the collisional processes do not conserve particle number and there- 
fore change chemical equilibrium, the only conserved quantity is the energy. Such is the 
case for massless particles interacting with a quartic couplings for example or higher order 
processes in a quartic theory with massive particles. Processes in which 3 — > 1 conserving 
energy and momentum can occur. The inverse process 1 — >• 3 occurs with far less probability 
since the high momentum modes are much less populated than the low momentum modes 
in the unstable bands. In this case only energy is conserved whereas the total number of 
particles is not conserved and in this case an estimate of the reheating temperature compares 
the energy density in the bath of produced particles to that of an ultrarelativistic gas in 
equilibrium at temperature Tr, 



\M. 



R\ 



l-l 



A; 

leading to the estimate 



T R (6.4) 



\M. 



R\ 



A fi 4 

Thus we can at least provide a bound for the reheating temperature 



(6.5) 



IMpI 

\M R \ <T r < (6.6) 

and a more quantitative estimate requires a deeper understanding of the collisional processes 
involved. 

Within the large N approximation, scattering terms will appear at order 1/N and beyond. 
The leading contribution 0(1/N) to collisional relaxation conserves particle number in the 
unbroken symmetry state because the product particles are massive. This can be seen from 
the fact that the self-energy to this order is given by the same chain of bubbles that gives 
the scattering amplitude but with two external legs contracted. All cut diagrams (that 
give the imaginary part) correspond to 2 — > 2 processes that conserve particle number 
because of kinematic reasons. Certainly at higher order in 1/N there will be processes that 
change chemical equilibrium, but for large N these are suppressed formally. This argument 
based on the leading collisional contribution in the 1/N expansion allows us to provide a 
further consistent estimate in the unbroken symmetry case (when the produced particles 
are massive). In this approximation and consistently with energy and particle number 
conservation we can assume that the final equilibrium temperature is of the order of the 
typical particle energy before thermalization. Recalling that the unstable band remains 
stable during the evolution with the peak shifting slightly in position we can estimate the 
typical energy per particle by the position of the peak in the distribution at qi, and use the 
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analytical estimate for the peak given in section 2. Restoring the units we then obtain the 
estimate 



T R « \M r \ qi 



\M R \ | - % 



(6.7) 



which displays the dependence on r] explicitly. Eq. ( |6.7[) is an improvement of the simple 
estimate (|6.3j). 

It must be noticed that the peak in the momentum distribution decreases with time for 
t > t re h (see fig.2f-2h). This drift follows from the non-linear interaction between the modes. 
For the case of fig. 2, one sees that T R reduces by approximately a factor 3 with respect to 
the value (|6.7|) . 

The large N model studied in this article is not a typical model used in inflationary 
cosmology, and since we want to make a quantitative statement for inflationary scenarios 
(within the approximation of only considering Minkowski space) we now study a model that 
incorporates other scalar fields coupled to the inflaton. 

The simplest model [f||| contains in addition to the inflaton a lighter scalar field a with 
a ga $ 2 coupling. That is, we consider the Lagrangian || , 



C 



(d 2 + m 2 + go) $ - 4 
2 v ; 4! 



$ 4 - -a {d 2 + ml) a 



Act 

4! 



a 



(6.8) 



We will consider again the preheating regime of weak couplings and early times such 
that we can neglect the back-reaction of the quantum fluctuations of the a field as well 
as the back-reaction of the quantum fluctuations of the inflation itself, focussing only on 
the parametric growth of the a fluctuations in the unbroken symmetry case. The mode 
equations for the a field take the form 



dt 2 



+ k 2 + m 2 a + g<j) 2 (t) 



V k (t) = 



(6.9) 



In dimensionless variables this equation becomes 



dr 2 



+ q z + 



m 



V k {r) = . 



(6.10) 



In the short time approximation we can replace t]{t) by the classical form ( |3.1| ). We then 

fits closec 

n(n + 1) 



find a Lame equation which admits closed form solutions for |18 

12 g 



A 



n 



1,2,3, 



(6.11) 



Although these are not generic values of the couplings, the solubility of the model and 
the possibility of analytic solution for these cases makes this study worthwhile. In the 
simplest case, (n = 1, 6g = A), there is only one forbidden band for q 2 > 0. It goes from 

<? 2 = 1 - (^f) 2 to q 2 = 1 - ) 2 + Vo/2. That is, 



forbidden band : m 2 — m 2 < k 2 < m 2 



m„ 



(6.12) 
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The Floquet index in this forbidden band takes the form 



F(q) = 2i K{k) Z(2K(k) v) ± tt , (6.13) 
where now v and q are related by the following equation 

= 1 _ (T£\ 2 + "A cn 2 (2K(k) v, k) . (6.14) 
V m J 2 

and < v < 1/2. 

The imaginary part of the Floquet index is now maximal at q\ = 1 — f^f 1 + and we 
can use this value to provide an estimate for the reheating temperature in this model in the 
same manner as for eq. fl6.7p , yielding the following estimate for the reheating temperature 



- reh 



1 m 2 - ml + - $q . (6.15) 



In the special case m a — m — \M R \, we recover eq.( |6.7|) , as expected. 

Again, the non- lineal field evolution for t > t re h decreases T re h. In the third reference 
under [||, we found for late times a T re h ten times smaller than the value fl6.15|) for g = 1.67T 2 . 

The estimates on the reheating temperature provided above should not be taken rigor- 
ously, but as an approximate guide. A consistent estimate of the reheating temperature and 
the thermalization time scales would in principle involve setting up a Boltzmann equation 
||. Under the assumption of a separation between the preheating and thermalization time 
scales one could try to use the distribution functions N q at the end of the preheating stage 
as input in the kinetic Boltzmann equation. However, we now argue that such a kinetic de- 
scription is not valid to study thermalization. A kinetic approach based on the Boltzmann 
equation, with binary collisions, for example, would begin by writing the rate equation for 
the distribution of particles 

N k (x\ 2 J d 3 k 1 d 3 k 2 d 3 k 3 <5 4 (£;i + k 2 + k 3 + k) [(1 + n k ){\ + n kl )n k2 n k3 - 

n k n kl (l + n fca )(l + n ks )]. (6.16) 

However this equation is only valid in the low density regime. In particular for the case under 
study, the occupation numbers for wavevectors in the unstable bands are non-perturbatively 
large oc 1/A and one would be erroneously led to conclude that thermalization occurs on the 
same time scale or faster than preheating. 

Clearly such a statement would be too premature. Without a separation of time scales 
the kinetic approach is unwarranted. The solution of the Boltzmann equation provides a 
partial resummation of the perturbative series which is valid whenever the time scales for 
relaxation is much longer than the microscopic time scales [27|, in this case that of particle 
production. 

In the case under study there is an expansion parameter, 1/N and clearly these scat- 
tering terms are sub leading in this formal limit, so that the separation of time scales is 
controlled. In the absence of such an expansion parameter, some resummation scheme must 
be invoked to correctly incorporate scattering. In particular when the symmetry is broken, 
the asymptotic excitations are Goldstone bosons, the medium is highly excited but with 
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very long-wavelength Goldstones and these have very small scattering cross sections. Such a 
resummation is also necessary in the large temperature limit of field theories in equilibrium. 
In this case the perturbative expansion of the scattering cross section involves powers of 
XT/m with m being the mass. A correct resummation of the (infrared) divergent terms 



leads to A(T) — > m/T in the large T/m limit ||28|| . In particular the 1/N corrections in the 



formal large N limit involve such a resummation, but in the non-equilibrium situation, the 
numerical implementation of this resummation remains a formidable problem. 



VII. CONCLUSIONS 

It is clear that preheating is both an extremely important process in a variety of settings, 
as well as one involving very delicate analysis. In particular, its non-perturbative nature 
renders any treatment that does not take into account effects such as the quantum back- 
reaction due to the produced particles, consistent conservation (or covariant conservation) of 
the relevant quantities and Ward identities, incapable of correctly describing the important 
physical phenomena during the preheating stage. 

In this work, we dealt with these issues by using the O(N) vector model in the large N 
limit. This allows for a controlled non-perturbative approximation scheme that conserves 
energy and the proper Ward identities, to study the non-equilibrium dynamics of scalar 
fields. Using this model we were able to perform a full analysis of the evolution of the zero 
mode as well as of the particle production during this evolution. 

Our results are rather striking. We were able to provide analytic results for the field 
evolution as well as the particle production and the equation of state for all these components 
in the weak coupling regime and for times for which the quantum fluctuations, which account 
for back-reaction effects, are small. What we found is that, in the unbroken symmetry 
situation, the field modes satisfy a Lame equation that corresponds to a Schrodinger equation 
with a two- zone potential. There are two allowed and two forbidden bands, which is decidedly 
unlike the Mathieu equation used in previous analysis |l|,|,12]. The difference between an 



equation with two forbidden bands and one with an infinite number is profound. We were 
also able to estimate analytically the time scale at which preheating would occur by asking 
when the quantum fluctuations as calculated in the absence of back-reaction would become 
comparable to the tree level terms in the equations of motion. The equations of state of both 
the zero mode and "pions" were calculated and were found to be describable as polytropes. 
These results were then confirmed by numerical integration of the equations, and we found 
that the analytic results were in great agreement with the numerical ones in their common 
domain of validity. 

When the O(N) symmetry is spontaneously broken, more subtle effects can arise, again 
in the weak coupling regime. If the zero mode starts off very near the origin, then the 
quantum back-reaction grows to be comparable to the tree level terms within one or at most 
a few oscillations even for very weak coupling. In this case the periodic approximation for 
the dynamics of the zero mode breaks down very early on and the full dynamics must be 
studied numerically. 

When numerical tools are brought to bear on this case we find some extremely interesting 
behavior. In particular, there are situations in which the zero mode starts near the origin 
(the initial value depends on the coupling) and then in one oscillation, comes back to almost 
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the same location. However, during the evolution it has produced 1/g particles. Given that 
the total energy is conserved, the puzzle is to find where the energy came from to produce the 
particles. We found the answer in a term in the energy density that has the interpretation 
of a "vacuuum energy" that becomes negative during the evolution of the zero mode and 
whose contribution to the equation of state is that of "vacuum". The energy given up by 
this term is the energy used to produce the particles. 

This example also allows us to study the possibility of symmetry restoration during 
preheating |p yi3l , p5| . While there have been arguments to the effect that the produced 
particles will contribute to the quantum fluctuations in such a way as to make the effective 
mass squared of the modes positive and thus restore the symmetry, we argued that they 
were unfounded. Whereas the effective squared mass oscillates, taking positive values during 
the early stages of the evolution, its asymptotic value is zero, compatible with Goldstone 
bosons as the asymptotic states. 

Furthermore, this says nothing about whether the symmetry is restored or not. This is 
signaled by the final value of the zero mode. In all the situations examined here, the zero 
mode is driven to a non-zero final value. At this late time, the "pions" become massless, i.e. 
they truly are the Goldstone modes required by Goldstone's theorem. 

The arguments presented in favor of symmetry restoration rely heavily on the effective 
potential. We have made the point of showing explicitly why such a concept is completely 
irrelevant for the non-equilibrium dynamics when profuse particle production occurs and 
the evolution occurs in a highly excited, out of equilibrium state. 

Finally, we dealt with the issue of how to use our results to calculate the reheating 
temperature due to preheating in an inflationary universe scenario. Since our results are 
particular to Minkowski space, we need to assume that preheating and thermalization occur 
on time scales shorter than the expansion time, i.e. H~ x . We also need to assume that 
there is a separation of time scale between preheating and thermalization. Under these 
assumptions we can estimate the reheating temperature as T re h oc \Mr\ in the case where 
the produced particles are massive and T reh oc |Mr|/A* in the massless case. 

We have made the important observation that due to the large number of long- wavelength 
particles in the forbidden bands, a kinetic or Boltzmann equation approach to thermalization 
is inconsistent here. A resummation akin to that of hard thermal loops, that consistently 
arises in the next order 'ml/N must be employed. In equilibrium such a resummation shows 
that the scattering cross section for soft modes is perturbatively small despite their large 
occupation numbers. 

There is a great deal left to explore. Clearly the first step from this work would be to 
generalize what we have done to include the expansion of the Universe. This should be 
done both in order to understand preheating in more detail, as well as to understand the 
evolution of the scalar field during the inflationary period. Further steps should certainly 
include trying to incorporate thermalization effects systematically within the 1/N expansion. 

As we were finishing the write-up of this work, we learned about complementary work by 
Cooper, Kluger, Habib and Mottola [[30] who studied similar issues in the broken symmetry 
phase with results that are consistent with those found by us. 
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APPENDIX A: SOME RESULTS ON ELLIPTIC FUNCTIONS 



Here we detail some of the derivations used in the main text for the analytic results. 



1. The Unbroken Symmetry Case 



We derive here the solutions of the mode equations 

+ q 2 + 1 + r,l - T] 2 sn 2 (ry/l + r%, fc) 



dr 2 



ip q {r) = 



(Al) 



where k 



■no 



It is convenient to express the Jacobi sine in terms of the Weierstrass V function through 



sn (rJl + rilk 



k 2 sn 2 [rJl + 1% + iK'(k),k 



1 



k 



k 2 (ei - e 3 ) 
•j e 2 - e 3 



V 



T^l + ri 2 + iK\k) 



e 3 



ei - e 3 



(A2) 
(A3) 



We assume that the discriminant A of the function V is here positive and that the roots 
ei, e 2 , e 3 of the cubic equation obeying e± + e 2 + e 3 = are ordered as 



e 3 < e 2 < < e x . 



(A4) 



In addition, without any loss of generality we choose e\ — e 3 = 1 + 7/ 2 . With such a choice 



we find the roots to be: 3ei = |r] 2 + 2, 3e 2 = — 1 and 3e 3 = —\r\ 2 — 1. 
Collecting all factors the mode equations (|A1|) become 

d 2 o 1 



dr 



2 • 1 



1 ' --2V(t + uj') 



if g (r) = 



(A5) 



where u/ 



■ K'(k) 
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functions. The solutions is 18 



Eg. ( |A5|) is a Lame equation that can be solved in closed form in terms of Weierstrass 

(A6) 



q{ ' a(r + u;')a(u' + z) 



where z is defined through 



Viz) 



(A7) 



<j{x) and ((x) are Weierstrass functions. We normalize the functions U q {r) according to 
eq.fj). 

Changing r by — r in eq . (|A6|) provides an independent solution of eq.( |A5|) . 
Eq. ( |A7| ) maps the real q 2 axis into the sides of the fundamental square in the z-plane 
MM. That is: 



z = p, < (3 < u , -oo < q 2 < -1 - 



2 ' 



4 



where uj 



K{k) 



z = u + ia , 0<a< > -1 - y < ? 2 < , 

2 

;z = o/ + /3, w>/5>0; 0<g 2 <^, 

2 

z = ia , u'/i >a>0; y<g 2 < +oo . 
. Notice that 2oj is the period in r of the potential in eqs.( |3.3|) and 



(A8) 



The mapping of the real g 2 -axis into the sides of the fundamental square in the z-plane is 
made explicit by writing the Weierstrass V(z) function in terms of Jacobi sn and cn functions 
LI23U ' 



v 2 + i 



VU' + iH) I -^- + (l+'f)sn 2 [a^l + V 2 ,k' 



P(^ + /3) = |cn 2 (/? v /l + % 2 ,M , 



(A9) 



- V(ia) 



-1-^1 + 
2 



v 2 + i 



sn 2 aJl + r/ 2 , k 



Here a and (3 are real variables. [Recall that < sn(w, k) < 1 for < u < K(k) and that 
sn(ii, k) 2 + cn(u, k) 2 = 1 .] 

Generically, a periodic potential possesses an infinite number of allowed and forbidden 
(Floquet) bands that alternate for < q 2 < +oo. A detailed study of the Floquet indices 
reveals that this is not the case for eq. (|Al|) . 
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Let us now compute the Floquet indices for the mode functions 



periodicity property of the Weierstrass cx-functions states that |20| , |23l , |24|1 : 

a(x + 2u) = -a(x) e 2{x+uj)r > 
where rj = ((uj). Using eqs.( [Al(i| ) and ( |3.4| ), we find from eq.( |A6| ) 

U q (r + 2u) = U q {r) e 2[zr >-^ z)] 

F{q) = 2i\u({z)-zrj\. 



. The quasi- 
(A10) 



(All) 



The quantity q belongs to an allowed band when 2[zt] — u((z)] is purely imaginary and 
q belongs to a forbidden band when this quantity has a nonzero real part. 

Using the properties of the Weierstrass functions [EDI], we find two allowed bands: 



and two forbidden bands 



z = ia and z = ia + u, 
with real a, < a < to' ji . 



z = P and z — (3 + u', 
with real (3, < (3 < lu . 



(A12) 



In terms of q 2 the allowed bands correspond to 



2 2 

1 - 5° < q 2 < and ^ < q 2 < +oo 



and the forbidden bands to 



oo<g 2 <-l-^ and < q 2 < ^ . 



(A13) 



(A14) 



(A15) 



The last forbidden band is for positive q 2 and hence does contribute to the fluctuation 
function S(r). S(r) will grow exponentially in time due to the presence of such unstable 
modes. 

2 

Let us investigate the Floquet indices for the forbidden band < q 2 < Setting 
z = j3 + to' we find 



2[zri-u({z)] 



2uj> 



VK 



Ml) 



(A16) 



where we used the relation between the Weierstrass ^-function and the Jacobi ^-functions 



C(z) 



(A17) 
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Here 77 = C(<^)- 

2 

In summary, we have on the forbidden band < q 2 < ^ , 



F(q) = ±tt + i 



±n + 2iK(k) Z(2K(k) v) 



(A18) 



where v = ^ , < v < 1/2 and Z(u) = E(u, k) — u is the Jacobi zeta function [IS 
Eq.( |3.12| ) gives its series expansion. We see that the solution U q (r) decreases with r. The 
other independent solution U q (—r) grows with r. 

Notice that U q (r) becomes just an antiperiodic function on the borders of this forbidden 
band, v — and v — 1/2. 

It is useful to write the solution U q {r) in terms of Jacobi ^-functions. Using eqs. 
( 15T7D and 



^ 2 ^ ^ Q / Z \ 



(A19) 



we found eq.( |3.10| ) after some calculation. 

For the allowed band ^= < q < 00, we find eq. (|3.21|) using eqs. (|A6|) , eq.( |A17| ) and 
eq. ([AT9D . 



2. Broken Symmetry Case 



We derive here the solutions of the mode equations fl4.5|), 



<p q (r) = 



(A20) 



It is convenient to express the Jacobi function dn in terms of the Weierstrass V function 

'u + K(k) +iK'(kY ~ 



through [EC 



1 



1 



dn 2 (-u, k) ei — e 2 



(A21) 



n 

Where now for convenience and without loss of generality we choose e\ — e% = 1 — f. Then, 
we find 3ei = 1 , 3e 2 = 1 — §r/o and 3e 3 = §r?o — 2 



Collecting all factors the mode equations (|A20| ) become 

tp,(T) = 



(A22) 



where u' = % K J^— . This equation is equivalent to eq. (IA5I) up to the sign in front of the ~ 

v 1 -^ — ' 

and a shift r — > r + uj. 



The solutions of eq. ( |A22j ) are given by 
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U g (r) 



er(r + uj + uj' + z)a{yj' + uj) 
cr(r + uj + uj') a{uj' + uj + z) 



(A23) 



2; and g 2 are now related by 



n*) = +3 - 1 2 



(A24) 



We continue to use the normalization ( |3.5| ). Changing r by — r in eq.( |A23| ) provides in 
general an independent solution of eq.( |A20|) . 

Eq. ( |A24|) maps the real q 2 axis into the four sides of the fundamental square in the z- 
plane. This mapping is better seen writing the Weierstrass V(z) function in terms of Jacobi 
sn and cn functions [EDI in the four cases. That is: 



z = p, < P <uj , -00 < g 2 < 0. 



1 



V((3) 



V 2 \ Cn 



/3\A 



sn' 



z = uj + ia , 0<a< uj' /i , < q 2 < y , 



1 



ni 



Vo 



z = uj' + (3,uj>(3>0; 'f<q 2 <l , } 
q 2 = \- V{J + (3) = 1 - I - (1 - V 2 )sn 2 I 

z = ia , u//i >a>0; l-y<g 2 < +00 

1 1 - 2a 

g 2 = --P(ia) 

sn 2 f a\/l - f,k' 



Vo 



where uj 



K(k) 
1 2 



(A25) 



. Notice that 2uj is the period in r of the potential in eqs. ( |A20| ) and 



Eq. (|3.4| ) holds for the mode functions ( |A23| ) too. Hence, there are two allowed bands: 

(A26) 



z = ia and z = ia + uj, 
with real a, < a < uj' /i . 



and two forbidden bands 



z = (3 and z = /3 + a/, 
with real (3, < (3 < uj . 



(A27) 



In terms of q 2 the allowed bands correspond to 
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0< g 2<I7o and !_M< 2< 



(A28) 



and the forbidden bands to 



— oo < q 2 < and 



< q < 1 - 



(A29) 



The last forbidden band is for positive q 2 and hence does contribute to the fluctuation 
function S(r). S(r) will grow exponentially in time due to the presence of such unstable 
modes. 

It is useful to write the solution U q (r) in terms of Jacobi -^-functions. For the forbidden 

2 2 

band % < q 2 < 1 — ^, we obtain eq. 



using eqs. ( |A23| ), ( |A17| ) and ( |A19| ), after some 
calculation. 

2 

For the allowed band 1 — ^ < q 2 < +oo , we find eq. ( [4.1 1| ) from eqs.( |A23| ), ( |A17| ) and 

2 

( |A19| ). The analogous expression ( |4 . 1 3| ) hold in the allowed band < q 2 < 
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FIGURES 

Figure 1: The ratio < p > /e for zero mode vs. \ R e /2\M R \ 4: for the unbroken 
symmetry case. 

Figure 2a: r](r) vs. r for the unbroken symmetry case with i] = 4, g = 1CT 12 . 

Figure 2b: (?£(t) for the same values of the parameters as in Fig. 2(a). The agreement 
with the analytic prediction is to within 5% for < r < 30. 

Figure 2c: gAf(r) for the same parameters as in fig. 2a. 

Figure 2d: gN q (r) vs. q for r = 40 for the same values of parameters as fig. 2a. 
Figure 2e:gN q (r) vs. q for r = 120 for the same values of parameters as fig. 2(a). 
Figure 2f: gN q (r) vs. q for r = 200 for the same values of parameters as fig. 2(a). 

Figure 2g: ( 2 |Mr| 4 ) ^( r ) ^ or ^ ne same values of the parameters as in Fig. 2(a). Asymp- 
totically the average over a period gives p^ ~ e/3. 

Figure 3: The ratio < p > /e for zero mode vs. \ R e /2\M R \ 4: for the broken symmetry 
case. 

Figure 4a: rj(r) vs. t for the broken symmetry case with rj = 10" 5 , g = 10~ 12 . 
Figure 4b: (?S(r) for the same values of the parameters as in fig. 4a. 
Figure 4c: gM{r) for the same parameters as in fig. 4a. 

Figure 4d: gN q (r) vs. q for r = 30 for the same values of parameters as Fig. 4a. 
Figure 4e: gN q (r) vs. q for r = 90 for the same values of parameters as Fig. 4a. 
Figure 4f: gN q (r) vs. q for r = 150 for the same values of parameters as Fig. 4a. 
Figure 4g: M 2 (t) vs. r for the same parameters as Fig. 4a. 
Figure 4h: e c /(r) vs. r for the same parameters as Fig. 4a. 
Figure 4i: en(t) vs. t for the same parameters as Fig. 4a. 
Figure 4j: Sc(t) vs. t for the same parameters as Fig. 4a. 

Figure 4k: ( 2 |m^ | 4 ) p{ T ) f° r the same values of the parameters as in Fig. 4a. Asymp- 
totically the average over a period gives p^ = e/3. 

Figure 5a: t)(t) vs. t for the broken symmetry case with r] = 10~ 2 , g = 10~ 5 . 
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Figure 5b: M 2 (t) vs. t for the same parameters as fig. 5(a). 
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TABLES 
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TABLE I. Quantum Fluctuations S(r) ~ e Br during the preheating period. 
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